Generation of Kerr Frequency Combs in Resonators with Normal GVD 
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We show via numerical simulation that Kerr frequency combs can be generated in a nonlinear 
resonator characterized with normal group velocity dispersion (GVD). We find the spectral shape 
of the comb and temporal envelope of the corresponding optical pulses formed in the resonator. 



Kerr optical frequency combs have recently attracted 
considerable attention because of their rich physics and 
their potential for multiple applications in science and 
technology 1]. Kerr combs are spontaneously generated 
from quantum fluctuations in solid state nonlinear opti- 
cal microresonators pumped with continuous wave light. 
Since the first demonstration of a comb observed in a 
fused silica microresonator the possibility of generat- 
ing combs in other types of resonators have been reported 
in numerous experiments. In accordance with experimen- 
tal observations, the observed combs preferably occur in 
resonators characterized with anomalous group velocity 
dispersion (GVD). This conclusion has also been con- 
firmed theoretically [3, 0] • Recent numerical simulations 
conclude that generation of Kerr combs is also possible 
and stable at all but zero GVD [5]. 

Generation of the Kerr comb results from modulation 
instability (MI) of the continuous wave light propagat- 
ing in a nonlinear host material of an optical resonator. 
The physical origin of the instability is the same as those 
observed in optical fibers, where MI is observed under 
conditions of anomalous GVD and positive cubic non- 
linearity. It is natural to expect that the resonant MI 
has properties similar to the non-resonant MI. However, 
it was noted that MI also takesplace in a fiber ring res- 
onator with a net normal GVD \G\ . The reason is that the 
resonant configuration introduces an additional degree of 
freedom in the system, namely frequency detuning of the 
pumping light from the eigenmode of the nonlinear res- 
onator [7]; this can shift the MI instability point to the 
region of normal GVD. 

There have been several experiments in which the 
comb formation was achieved in microresonators with 
normal GVD. For instance, a comb was observed in a 
CaF2 whispering gallery mode (WGM) resonator with a 
diameter of 0.78 cm, pumped with 1320 nm light, and 
characterized with /32 — 7.7 ps^/km normal GVD [g]. A 
comb was also generated in a CaF2 WGM resonator fea- 
turing a diameter of 0.255 cm, pumped with 1550 nm 
light, and having (32 — 2.3 ps^/km Q. 

The experimental results were explained theoretically. 
It was shown that the resonant hyper-parametric oscilla- 
tion can occur at any GVD (the term "hyper-parametric 
oscillation" was introduced to describe the underdevel- 
oped Kerr frequency combs, containing only two sym- 
metric sidebands with respect to the carrier) [13, . 

The excitation dynamics represents the basic difference 
in Kerr combs generated in resonators with normal and 
anomalous GVD. Kerr combs can have both hard and 



soft oscillation onset in the case of anomalous GVD, but 
only the hard oscillation onset takes place in the case of 
normal GVD, as noted in [ll|. Nevertheless, a more rig- 
orous study of the subject is still required. The theory in 
prj was developed for a hyper-parametric oscillator only, 
and does not cover generation of Kerr frequency combs. 
The numerical simulations [Ij do not describe hard exci- 
tation of the comb and do not allow finding any steady 
state solution for the comb, since it is assumed that the 
generated sidebands are always smaller compared to the 
optical pump. 

In this Letter we present results of our numerical sim- 
ulations of Kerr frequency combs produced in a non- 
linear microresonator with normal GVD. For the sake 
of specificity, we present a steady state solution for the 
comb generated in a GaF2 WGM resonator pumped at 
1554 nm. The resonator is characterized with 667 /im 
diameter, 100 GHz FSR, no = 1.43 refractive index, 
and (32 — 0.05 5 ps^/km GVD, corresponding to 110 kHz 
frequency difference between adjacent FSRs {v+ + — 
2vf) ~ —110 kHz, where vq is the linear frequency of the 
pumping light, v± are the frequencies of the optical side- 
bands, — ~ vq — v- ~ 100 GHz). The resonator host 
material has cubic nonlinearity 712 = 1.9 x 10~^^ cm^/W, 
and mode volume V = 1.3 x 10^^ cm^^. We assume that 
the full width at the half maximum of the overloaded 
WGMs is 270 = 220 kHz, so that v+ + v_- 2iyo = -70- 

To find the GVD we used CaF2 Sellmeier equation 
[12] and an asymptotic expression describing the spec- 
trum of a dielectric spherical resonator 13]. The over- 
all dispersion /?2(t'o) was calculated using the formula 
iy+ + L'-- 2i/o = -7rc/32(j^o)(j^+ - I'- f /2no{L'o) M- The 
selected value of GVD {/32) seems to be small, since it is 
much smaller compared with the dispersion of a conven- 
tional optical fiber, however the significance of GVD is 
determined by the ratio {v^ + iy- — 2i/o|/7o in a resonator 
p^ . The GVD is considered small if this parameter is 
much less than unity. Since the difference + p_ — 2pq 
depends on the FSR, larger resonators have smaller rel- 
ative GVD than smaller resonators. Resonators with 
higher Q-factor have larger relative dispersion compared 
with lower Q resonators. 

To study the behavior of the comb numerically it 
is convenient to introduce a coupling constant [To| 
g = hujlcn2 / (Vnl) = 3.3 x 10"^ s"\ where ujq = 
2-kvq. Selecting the pump power P — 36.7 /iW, 
we also introduce the dimensionless pumping constant 
/ = (Fo/27r7o)(5/27r7o)^/^ = 2.244, where Fo = 
{ATT^fjP / {Tlloq))'^/'^ describes the amplitude of the contin- 
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uous wave external pump. The selection of the above 
listed parameters is an arbitrary one, simply serving the 
goal of comparison of the results obtained in our numer- 
ical simulations with parameters of a realizable system. 

In our model we numerically study the nonlinear in- 
teraction of 21 optical modes. The selected number of 
modes is limited by the available computational capac- 
ity. We expect that all the modes are identical and com- 
pletely overlapping in space. The external continuous 
wave pump is applied to the central mode of the mode 
group, so the simulated Kerr comb is expected to have 
ten red- and ten blue-detuned harmonics with respect to 
the frequency of the pumped mode. We take into ac- 
count only the second order frequency dispersion that is 
recalculated to the frequency of the modes 

To write the nonlinear equations we introduce an in- 
teraction Hamiltonian V — —g{e^)^e^ /2, where e = 
S^=i ^j- The equations of motion are [sl. [lOj 

= -(70 + i^j)aj + ^[V, a^] + Foe-*"*(5iij, (1) 

where Su j is the Kronecker's delta. The set of nonlin- 
ear forced equations, ([l}, is solved numerically without 
any further assumptions. The evaluation was interrupted 
when the solution reached its steady state (70*3 ~ 30). 
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FIG. 1: Normalized slow amplitude (|^ii|(t;/27r7o)^''^) and 
phase of the intracavity field of the pumped mode ver- 
sus normalized detuning ((cj — a)ii)/27r7o). The region of 
the maximum field accumulation is shifted to smaller fre- 
quencies due to the self-phase modulation effect (the effec- 
tive index of refraction increases with intensity of light cir- 
culating in the mode). The pump power is fixed, / — 
(Fo/27r7o)(ff/27r7o)^''^ = 2.244. The stable solutions of the 
set {TJ that do not reveal the Kerr comb generation are shown 
by blue lines. The stable solution that reveals the comb gen- 
eration is shown by red line. The solution is localized in both 
amplitude and frequency spaces. The inset shows the region 
of parameters where the comb is generated. 

The results of simulations are summarized in Figs. ([T])- 
(|3]). To present the results we introduce slow amplitudes 
of the field in the modes dj — Aj exjp^—iojjt). The slow 
amplitudes are further normalized as \Aj\{g/2Tr^oy^^ . 
The frequency detuning is also normalized to the half 
width at half maximum of the resonator mode, ((w — 
a;ii)/27r7o). As a rule, the solution followed stable 
branches (shown by solid blue line) and then jumped to 



the attractor corresponding to comb generation (shown 
by the red line). The jump was observed only for the 
case of nonzero initial conditions. Naturally, the solution 
does not converge to unstable solutions. 

We find that a Kerr comb can be generated in the 
resonator in the case of normal GVD. For the resonator 
described above we obtain the power and the frequency 
of the pumping light that result in comb excitation. The 
frequency comb is dynamically stable, and the genera- 
tion onset is hard in its nature, i.e. there is a discontinu- 
ous jump in parameters of the comb harmonics when the 
pump power exceeds a certain threshold. 
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FIG. 2: Normalized slow amplitude {\Aii\{g/2TT-yoy^^) and 
phase of the intracavity field of the pumped mode ver- 
sus the normalized amplitude of the external pump (/ — 
{Fo/2TT'yo){g/2n'yo)^^^)- The frequency detuning is fixed 
(id — ci;ii)/27r7o = —4. Stable solutions of the set fl]) that 
do not reveal Kerr comb generation are shown by blue lines. 
The stable solution that reveals comb generation is shown by 
the red line. The solution is localized. The inset shows the 
region of parameters where the comb is generated. 

The behavior of the field of the mode pumped optically 
is illustrated by Figs. ([1]) and ([2]). If the comb sidebands 
are small the stable oscillation occurs in the region that is 
unstable. This stability region was not found in previous 
research 5] since only solutions with adiabatically stable 
amplitude of the pump mode were analyzed there. The 
corresponding behavior of amplitudes of the first two op- 
tical sidebands (^10 and A12) is shown in Fig. ([3]). The 
stability region where the sidebands are generated is lo- 
calized with respect to the frequency and power of the 
pumping light. 

We found the time dependence of the overall field am- 
plitude in the resonator (Fig. |3]) by selecting specific val- 
ues of power and frequency of the pump light. Appar- 
ently, the modes of the frequency comb create a single 
optical pulse traveling in the resonator. The shape of 
the pulse is rather complex and more effort is required 
to derive an analytical expression for its description. We 
verified that if the power of the external pump is changed 
a similar pulse shape, but with different amplitude and 
duration, is generated in the resonator . 

The initial conditions and parameters of equations are 
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FIG. 3: Normalized slow amplitude (|j4io|((//27r7o)^''^ and 
|j4i2|((7/27r7o)^''^ are the same) of the first sidebands of the 
pump mode for the same parameters of the system as used in 
Figs. (O and jH). 



varied to verify the stability of the solution. We also re- 
duced and increased the number of modes in our code 
and confirmed that the spectrum does not change signif- 
icantly. It is worth noting that the method is suitable for 
the simulation of a resonator with arbitrary GVD. The 
limited number of modes used in the simulation imposes 
a practical restriction on the values of the pump power 
and GVD values on the outcome of the simulations. We 
intentionally selected a large value of the relative GVD 
to suppress this undesirable effect. 
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FIG. 4: Normahzed slow amplitude {{g\e^ e\/2-K'yoy^'^) of the 
intracavity field versus time, and the corresponding optical 
comb found for a fixed pump power / = 2.244 and fre- 
quency detuning of the pump light from the selected mode 
(cj - a;ii)/27r7o = -4. The 100 GHz FSR of the resonator is 
expected. 



To conclude, we have studied Kerr frequency comb 
generation in microresonators made out of a transparent 
optical material possessing positive cubic nonlinearity, as 
well as normal group velocity dispersion. We find that 
there exists a stable attractor for oscillation that results 
in comb generation. Stable Kerr frequency combs can 
be produced in such resonators if the proper power and 
frequency for the continuous wave light pumping a mode 
of the resonator is selected. 
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